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ABSTRACT:  Direct  dynamics  simulations  and  selected  ion  flow  tube  (SIFT)  experiments 
were  performed  to  study  the  kinetics  and  dynamics  of  the  OH~  +  CH3I  reaction  versus 
temperature.  This  work  complements  previous  direct  dynamics  simulation  and  molecular 
beam  ion  imaging  experiments  of  this  reaction  versus  reaction  collision  energy  (Xie  et  al. }. 

Phys.  Chem.  A  2013,  117,  7162).  Tlie  simulations  and  experiments  are  in  quite  good 
^eement  Both  identify  the  Sn2,  0H~  +  CH3I  — >  CH3OH  +  I",  and  proton  transfer, 

OH~  +  CH3I  -*  CH2r  +  H2O,  reactions  as  having  nearly  equal  importance.  In  the 
experiments,  the  Sn2  pathway  constitutes  0.64  ±  0.05,  0.56  ±  0.05,  0.51  ±  0.05,  and  0.46 
±  0.05  of  the  total  reaction  at  210,  300,  400,  and  500  K,  respectively.  For  the  simulations 
this  fraction  is  0.56  ±  0.06,  0.55  ±  0.04,  and  0.50  ±  0.05  at  300,  400,  and  500  K, 
respectively.  The  experimental  total  reaction  rate  constant  is  (2.3  ±  0.6)  X  lO”’,  (1.7  ± 

0.4)  X  lO"*,  (1.9  ±  0.5)  X  lO"’,  and  (I.8  ±  0.5)  X  lO"*  cm* s"*  at  210,  300,  400,  and  500 
K,  respectively,  which  is  approximately  25%  smaller  than  the  coflision  capture  value.  The 
simulation  values  for  this  rate  constant  are  (1.7  ±  0.2)  X  lO*®,  (I.8  ±  O.l)  X  lO”’,  and  (I.6 
±  0.1)  X  10“*  cm*s~*  at  300,  400,  and  500  K.  From  the  simulations,  direct  rebound  and  stripping  mechanisms  as  well  as  multiple 
indirect  mechanisms  are  identified  as  the  atomic  level  reaction  mechanisms  for  both  the  Sn2  and  proton  transfer  pathways.  For 
die  5)^2  reaction  the  direct  and  indirect  mechanisms  have  nearly  equal  probabilities;  the  direct  mechanisms  are  slightly  more 
probable,  and  direct  rebound  is  more  important  than  direct  striping  For  the  proton  transfer  pathway  the  indirect  mechanisms 
are  more  important  than  the  direct  mechanisms,  and  striping  is  significantly  more  important  than  rebound  for  the  latter. 
Calculations  were  performed  with  the  OH~  quantum  number  /  equal  to  0,  3,  and  6  to  investigate  the  efiect  of  OH“  rotational 
excitation  on  the  OH~  +  CH3I  reaction  dynamics.  The  overall  reaction  probability  and  the  probabilities  for  the  S)^2  and  proton 
transfer  pathways  have  little  dependence  on  /.  Possible  effects  on  the  atomistic  mechanisms  were  invest^ated  for  the  S|^2 
pathway  and  the  probability  of  the  direct  rebound  mechanism  increased  with  /.  However,  the  other  atomistic  mechanisms  were 
not  appreciably  affected  by  /. 
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I.  INTRODUCTION 

Bimolecular  reactions  of  an  anion  X”  with  a  methyl  halide 
CH3Y,  e.g,  the  8^2  reaction  X~  +  CH3Y  CH3X  +  Y~,  are 
widely  studied  and  are  important  ion— molecule  reactions, 
particularly  in  organic  chemistry.*  The  dynamics  of  these 
reactions  have  been  extensively  studied  bodi  experimentally*”^* 
and  computationally.*”^  During  the  past  few  years  it  has 
become  possible  to  directly  observe  the  dynamics  of  ion- 
molecule  reactions  by  kinematically  complete  imaging  of 
reactive  scattering,  which  has  been  applied  to  charge  transfer,®’’ 
proton  transfer,  and  8^2  bimolecular  nucleophilic  sub 
stitution  reactions.*’”**  These  experiments,  along  with  direct 
dynamics  simulations, *'*’’*  *'**”**  have  provided  an  atomistic 
understanding  of  the  dynamics  of  gas  phase  X”  +  CH3Y 
reactions. 

In  recent  work,*®  ion  imaging  and  direct  dynamics 
simulations  were  used  to  study  the  atomistic  mechanisms  for 
the  OH”  +  CH3I  reaction  at  collision  energies  of  0.05,  0.5,  1.0, 
1.5,  and  2.0  eV.  Direct  comparisons  between  the  experiments 
and  simulations  were  made  at  0.5,  1.0,  and  2.0  eV,  and  quite 


good  agreement  was  found.  However,  no  comparison  between 
experiment  and  simulation  was  possible  for  the  lowest  collision 
energy  studied  in  the  simulations  of  0.05  eV,  which  is 
representative  of  a  temperature  of  ~400  K.  We  have  therefore 
undertaken  a  temperature  dependent  study  of  the  reaction  of 
OH”  with  CH3I  from  210  to  500  K,  by  both  e^qjeriment  and 
direct  dynamics  simulations. 

The  research  reported  here  addresses  atomic  level  mecha 
nisms  for  the  OH”  +  CH3I  reaction.  The  traditional  pathway 
for  X”  +  CH3Y  reactions  is  the  5n2  reaction,  for  which  X” 
CH3Y  and  XCH3  Y”  pre  and  postreaction  complexes  are 
formed  and  separated  by  a  central  barrier.  Both  simulation  and 
experiments  have  shown  that  X”  +  CH3Y  reactions  are 
considerably  more  rich  in  detail  with  additional  pathways  and 
numerous  atomistic  mechanisms,  as  summarized  in  Table  1. 
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Table  1.  X  +  CHjY  S|^2  and  Proton  Transfer  Reaction  Mechanisms 

Direct  RebounJ:  X~  attacks  the  backsides  of  CH3Y,  directly  replacing  Y~  for  Sj42  (or  abstracting  an  H-atom  for  proton  transfer)  reaction,  and  scattering  in  the 
backward  direction,  with  inversion  for  die  reaction.  Tends  to  be  a  small  impact  parameter  event;  refe  11  and  18. 

Stripping:  X~  approaches  the  side  of  CH3Y,  stripping  off  CHj  for  5)^2  (or  an  H-atom  for  proton  transfer)  reaction,  and  scattering  in  the  forward  direction, 
iiith  CH3  inversion  for  the  $1^2  reaction.  Tends  to  be  a  large  impact  parameter  event;  refe  11  and  18 

Front  Side  Attack:  X~  attacks  the  front  side  of  CHjY  and  directly  replaces  Y~  widiout  CH3  inversion;  ref  18. 

Indirect  Ion  Dipole  Complex:  X~ — CH3Y  and  XCH3 — Y~  pre-  and  postreaction  complexes;  refs  2  and  5. 

Hydrogen-Bonded  Complex:  X~ — HCHjY  prereaction  and  CH3X— Y“  (X  =  OH“)  postreaction  complex;  refe  11  and  18 
Proton-Tranter  Complex:  HX — CHjY^  (X  =  OH“);  this  work.  [CH3 — Y — X]“  Complex-,  refe  18  and  19. 

Roundabout  XT  collides  with  CH3Y  and  the  CH3-group  rotates  aroimd  Y  one  or  more  times  before  Sn2  substitution  or  proton  transfer  occurs;  ref  10. 
Proton  Exchange:  XT  =  OH“  mechanism  for  which  there  is  an  exchange  between  the  H-atom  of  OH“  and  a  H-atom  of  Ol3Y;  ref  18. 

Barrier  Recrossing:  nonstatistical  recrossing  of  the  5^2  [X— CHj-Y]“  central  barrier;  refe  14,  20,  and  21. 


II.  COMPUTATIONAL  PROCEDURE 

The  computational  procedure  used  for  the  simulations  was 
described  in  detail  previously*®  and  is  briefly  summarized  here. 

A.  Electronic  Structure  Calculations.  The  NWChem 
computer  program was  used  for  the  electronic  structure 
calculations.  Density  functional  theory  (DFT)  calculations,^ 
with  the  B97  1  functional^  and  the  ECP/d  basis  seL  give 
reaction  energies  in  good  agreement  with  experiment.*®  For  the 
ECP/d  basis  set,^  Dtmning  and  Woon’s  av^  cc  pVDZ  basis 
^^27,28  used  for  the  H,  C,  and  O  atoms.  For  iodine,  the 
Wadt  and  Hay  efliective  core  potential  (ECP)^^  was  used  for 
the  core  electrons  and  a  3s,  3p  basis  set  for  the  valence 
electrons,  which  was  augmented  by  a  d  polarization  function 
with  a  0.262  eiqionent  and  s,  p,  and  d  difluse  fimctions  with 
exponents  of  0.034,  0.039,  and  0.0873,  respectively. 

B.  Direct  Dynamics  Simulations.  Tlie  direct  dynamics 
simulations®**  for  the  OH“  +  CH3I  reaction  were  performed  at 
the  B97  1 /ECP/d  level  of  theory.  The  trajectories  for  the 
simulations  were  calculated  using  the  VENUS  chemical 
dynamics  computer  program®*"®®  interfaced  with  the  NWChem 
electronic  structure  computer  program.®®"®®  To  directly 
compare  with  the  experiments,  the  simulations  were  performed 
for  reactant  temperatures  T  of  300,  400,  and  500  K.  Initial 
conditions  for  the  trajectories  were  selected  for  collision 
energies  (Le.,  reagent  relative  translational  energy)  of  0.90,  1.05, 
and  1.20  kcal  mol”*,  equivalent  to  3RT/2  for  the  temperatures 
of  3(X),  400,  and  500  K.  Quasiclassical  sampling  for  temperature 
T  was  used  to  select  initial  conditions  for  the  CH3I  vibrational 
levels,  while  the  thermal  rotational  energy  of  CH3I  was  treated 
classically.  Thermal  sampling  was  used  to  sample  the  OH” 
rotational  and  vibrational  quantum  numbers.  The  OH” 
rotational  quantum  number  /  was  selected  according  to  the 
distribution  P(J)  =  (2/  -f  1)  exp[-/(/  -b  l)h^/2IkiT]/Qjof  The 
vibrational  quantum  number  n  was  selected  from  the 
distribution  P{n)  =  [exp(  nhv/k^T)]/Q^ 

In  previous  simulations,*®  the  OH”  -b  CH3I  reaction  was 
studied  at  a  collision  energy  of  0.05  eV  for  fixed  OH”  rotational 
and  vibrational  quantum  numbers  /  =  3  and  n  =  0.  The  CH3I 
rotational  and  vibrational  temperatures  were  T,  =  130  K  and  T, 
=  330  K  To  study  how  OH”  rotational  energy  afiects  the 
reaction  dynamics,  simulations  were  performed  with  these 
initial  conditions  but  with  /  equal  to  0  and  6,  to  compare  with 
the  previous  /  =  3  simulation  results. 

Semiclassical  quantization  of  the  action  integral®®  for  the 
potential  [V(r)  +  /(/  +  l)^®/(2//r®)]  was  performed  to  select 
initial  conditions  for  the  OH”  n,  J  vibrational  and  rotational 
energy  levels.  Algorithms  for  sampling  the  OH”  +  CH3I  initial 
conditions  are  standard  options  in  VENUS  and  have  been 
described  previously.®^ 
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A  sixth  order  symplectic  algorithm®®"®^  with  a  25  fe  time  step 
was  selected  to  integrate  the  direct  dynamics  trajectories.  The 
total  integration  time  for  the  trajectories  was  2.5  ps.  Reactive 
trajectories  were  identified  by  viewing  animations  and  thus 
determining  their  atomic  level  mechanisms.  The  properties 
determined  from  the  simulations  are  as  follows:  reaction 
pathways  and  their  atomic  level  mechanisms  and  reaction  cross 
sections  for  the  difierent  pathways  and  their  individual 
atomistic  mechanisms.  The  above  properties  compared  with 
tfie  experiments  of  \Tggiano  and  co  workers  are  the  relative 
product  ion  yields  and  the  reaction  cross  sections  and  rate 
constants  for  the  difierent  reaction  pathways. 

III.  EXPERIMENTAL  METHOD 

The  experiments  were  performed  using  the  Air  Force  Research 
Laboratory  variable  temperature  selected  ion  flow  tube  (SIFT). 
The  instrument  has  been  well  described  previously,®®  and  only 
details  important  to  the  present  measurements  are  given.  OH” 
was  made  in  a  higji  pressure  ion  source  through  electron  impact 
on  H2O.  The  OH”  was  mass  selected  using  a  quadmpole  mass 
filter,  injected  into  die  flow  tube  through  a  Venturi  inlet,  and 
carried  downstream  by  a  helium  bufier.  After  typically  3  ms  and 
qjproximately  10'*  collisions  of  the  OH”  with  helium,  CH3I  was 
added  through  a  finger  inlet.  After  a  typically  4  ms  reaction 
time,  the  core  of  flie  flow  was  sampled  flirough  a  small  orifice  in 
a  rounded  nose  cone,  while  the  bulk  of  the  gas  was  pumped  by 
an  oil  free  Roots  type  pump.  Sampled  ions  were  injected  into  a 
quadmpole  mass  spectrometer  and  detected  using  a  discrete 
dynode  particle  multiplier,  operated  in  ion  counting  mode. 

An  flows  were  measured  using  MKS  mass  flow  controners. 
The  ion  velocity  was  calculated  from  the  helium  flow  velocity 
and  a  previously  measured  ion  to  helium  velocity  ratio  of 
~1.6— 1.7.  Coupled  with  the  known  reaction  length  of  59  cm, 
this  yields  the  reaction  time.  Rate  constants  were  obtained  by 
monitoring  the  primary  decay  of  OH”  as  a  fimction  of  added 
neutral  concentration  in  the  flow  tube  as  shown  in  Figure  la  for 
a  500  K  data  set.  Decays  were  exponential  (linear  on  die 
semilog  plots  shown)  over  at  least  1  order  of  magnitude. 
Absolute  imcertainties  in  the  rate  constants  are  estimated  to  be 
±25%.®®  Branching  ratios  were  determined  by  extrapolating 
plots  of  percent  fraction  of  each  product  versus  neutral 
concentration  to  zero  neutral  flow  as  shown  in  F^re  lb.  When 
the  CH3I  concentration  was  increased,  the  I”  fraction  increased 
at  the  expense  of  CH2r,  showing  that  CH2I”  reacted  with 
CH3I  to  form  I”.  At  low  OH”  depletion  the  total  count  rate 
decreased,  while  at  h^her  levels  of  depletion  the  total  counts 
were  constant  within  about  10%,  even  though  CH2I”  was  being 
converted  into  I”.  This  behavior  indicates  that  there  is  some 
detection  discrimination  against  OH”,  while  the  discrimination 
between  I”  and  CH2I”  was  small.  Therefore,  we  estimate  the 
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Figure  1.  (a)  Ion  counts  as  a  fiinction  of  CH3I  concentration  at  500  K 
for  the  reaction  of  OH”  with  CH3I,  with  the  best  fit  line  to  the  primary 
ion  decay  shown,  (b)  Fraction  of  products  as  a  function  of  CH3I 
concentration;  branching  fi-action  values  are  obtained  by  extrapolating 
to  a  CH3I  concentration  of  zero. 


product  fractions  to  be  accurate  to  about  0.05  absolute  (10% 
discrimination  for  equal  products)  and  0.03  relative. 

IV.  RESULTS  AND  DISCUSSION 

A.  Pathways,  Anion  Product  Ratio,  and  Rate  Con¬ 
stants.  7.  Simulation  Results.  For  300,  400,  and  500  K,  for 
which  the  cottision  energy  is  3RT/2,  a  respective  total  of  311, 
385,  and  301  trajectories  were  calculated,  and  the  impact 
parameter  was  randomly  selected  from  0  A  to  the  maximum 
impact  parameter  of  9.5,  9.2,  and  9.0  A,  respectively.  The 
was  identified  by  the  impact  parameter  at  which  no 
reaction  occurred  out  of  100  trajectories.  The  cross  section  for 
a  reaction  pathway  cr,  was  obtained  by  multiplying  xb^ma 
reaction  probability  P,  for  the  pathway.  The  reaction  rate 
constant  fc  is  determined  by  multiplying  the  overall  total 
reaction  cross  section  c,  by  the  averse  relative  translational 
velocity 


The  following  three  pathways  were  observed  for  each  of  die 
three  temperatures,  and  their  reaction  probabilities  are  shown 
in  Table  2. 


OH”  +  CHjI  -  CH3OH  +  r 

(1) 

OH”  +  CHjI  CHjF  +  H2O 

(2) 

OH”  -I-  CH3I  -*  [CH,-I-OH]” 

(3) 

The  overall  reaction  probability  ranges 

from  0.42  to  0.44, 

varying  very  little  for  the  three  temperatures.  The  Sn2  pathway 
1  has  a  higgler  probability  than  that  for  proton  transfer  2  at  a  T 
of  300  and  400  K,  but  has  a  lower  probability  at  500  K.  The 
probability  of  the  5^2  pathway  is  0.23,  0.23,  and  0.20,  for  300, 
400,  and  500  K,  respectively.  The  corresponding  probabilities 
for  the  proton  transfer  pathway  are  0.18,  0.20,  and  0.21.  This  is 
the  same  trend  as  found  in  the  previous  study;'®  Le.,  a  high 
collision  energy  gives  more  CHjI”  +  HjO,  while  a  low  collision 
energy  gives  more  CH3OH  +  I”.  As  discussed  in  the  previous 
study,'®  the  [CH3  1  OH]”  ion  in  reaction  3  is  an  intermediate 
which  has  sufiScient  energy  to  dissociate  if  the  trajectories  were 
run  for  a  longer  period  of  time.  It  can  either  dissociate  back  to 
reactants  CH3I  +  OH”,  or  react  to  follow  the  5^2  or  proton 
transfer  pathways. 

2.  Comparison  between  Experiment  and  Simulation.  The 
diree  anions  I”,  CHjI”,  and  [CH3  1  OH]”  are  formed  in  die 
simulations  and  their  product  ratios,  given  in  Table  3,  are 
obtained  directly  from  the  pathway  probabilities.  Only  I”  and 
CH2I”  are  detected  in  the  experiments,  and  their  ratio  is  also 
given  in  Table  3.  Experimentally  the  I”  yield  decreases  from 
0.64  to  0.46  as  the  temperature  is  increased  fiom  210  to  500  K. 
This  trend  is  also  identified  fiom  the  simulation. 

The  [CH3  I  OH]”  anion  observed  in  the  simulations  at  the 
2.5  ps  conclusion  of  the  tr^ectory  integratiotis  will  dissociate  at 
longer  times.  An  incomplete  study  indicates  that  it  dissociates 
preferentially  to  the  reactants  HO”  +  CH3L  If  it  is  assumed  that 
all  of  the  dissociation  is  to  HO”  +  CH3I,  not  affecting  die 
relative  CH2I”  to  I”  yield,  the  fraction  of  the  I”  anion  in  die 
simulations  becomes  0.56,  0.55,  and  0.50  at  300,  400,  and  500 
K,  respectively. 

The  total  reaction  rate  constant  is  another  important 
property  that  can  be  compared  between  the  experiments  and 
simulations.  As  shown  in  Table  4,  the  experimental  and 
^Illation  rate  constants  are  in  quite  good  agreement.  Also 
given  is  the  collision  capture  rate  constant^®  for  each 
temperature,  using  the  experimental  data  for  the  OH”  + 
CH3I  km— dipole  interaction.  The  experimental  rate  constant  is 
^proximatety  25%  smaller  than  this  capture  value. 

The  simulation  rate  constants  are  based  on  the  formation  of 
the  CH2I”,  I”,  and  [CH3  1  OH]”  anions.  The  latter  species  is 
unstable  and,  as  discussed  above,  is  expected  to  preferentially 
dissociate  to  the  OH”  +  CH3I  reactants.  If  it  is  removed  from 


Table  2.  Temperature  Dependence  of  the  Simulation  Probabilities  for  the  OH  +  CH3I  Reaction  Pathways” 


probability 


temperature  (K) 

total  reaction 

CH3OH  +  r 

CHir  +  HjO 

[CH,-I-OH]- 

300 

0.42  ±  0.03 

0.23  ±  0.02 

0.18  ±  0.02 

001  ±  0.01 

400 

0.44  ±  0.03 

0J3  ±  a02 

0  JO  ±  0.02 

0.01  ±  0.01 

500 

0.43  ±  0.03 

0.20  ±  0.02 

OJl  ±  0.02 

002  ±  0.01 

“The  trajectories  were  terminated  at  an  integration  time  of  2.S  ps.  At  longer  times  all  of  the  [CH3  I  OH]”  anion  will  dissociate.  An  incomplete 
study  indicates  it  will  preferentially  dissociate  to  the  CH3I  +  OH”  reactants. 
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Table  3.  Comparison  of  the  Experiment  and  Simulation  Temperature  Dependence  of  the  Anion  Product  Ratio'’’*’ 


temperature  (K) 

CHjI- 

[CHj-l-OH]- 

r 

r  (exp) 

CHjI-  (eip) 

210 

0.64 

036 

300 

0.43  ±  0.05 

0.03  ±  0.01 

0.54  ±  0.06 

0.56 

0.44 

400 

0.45  ±  ao4 

0.01  ±  0.01 

a54  ±  0.04 

031 

a49 

500 

0.48  ±  0.06 

0.05  ±  0.01 

a47  ±  0.05 

0.46 

034 

“See  footnote  a  in  Table  2  and  IVAl  Simulation  Results  concerning  the  [CH3  I  OH]”  yield.  *Txperimental  branching  error  estimated  to  be  0.03 
relative  and  0.05  absolute. 


Table  4.  Comparison  of  the  Experiment  and  Simulation 
Temperature  Dependence  of  the  OH”  +  CHjl  Rate 
Constant"  and  the  Calculated  Collision  Rate  Constant 


rate  constant  (cmVs) 


temperature  (K) 

simulation 

exp"’ 

calc  coU.‘^ 

210 

23  X  10^ 

3.0  X  10^ 

300 

1.8  ±  0.1  xlO”’ 

1.7  X  10^ 

2.7  X  10^ 

400 

1.8  ±  0.1  xlO”* 

1.9  X  lO”* 

2.5  X  10^ 

500 

1.7  ±  0.1  xlO”* 

1.8  X  10^ 

2.3  X  10^ 

“If  the  [CH3  I  OH]”  anion  is  not  included  in  the  product  yield  (see 
footnote  a  in  Table  2),  the  simulation  rate  constant  becomes  1.7  X 
10”’  1.8  X  lO”’,  and  1.6  X  10”*  cmV  at  300,  400,  and  500  K, 
respectively.  ^’Errors  in  the  eiqjerimental  rate  constants  are  estimated 
to  be  ±25%.^^  The  collision  rate  constants  are  calculated  according  to 
tef  38. 


the  product  yield,  the  simulation  rate  constant  becomes  1.7  X 
lO”*,  1.8  X  10”’,  and  1.6  X  10”*  cm^  s”‘  at  300,  400,  and  500  K, 
respectively. 

B.  Potential  Energy  Surface  for  the  Proton-Transfer 
Pathway.  To  assist  in  understanding  the  atomistic  mechanisms 
for  the  OH”  +  CH3I  -»■  H2O  +  CH2I”  proton  transfer  reaction, 
detailed  electronic  structure  calculations  at  the  B97  1/ECP/d 
level  of  theory  were  performed  to  investigate  the  reaction’s 
stationary  points,  i.e.,  potential  minima  and  saddle  points. 
Figure  2  presents  their  structures  and  relative  energies  without 
zero  point  energy.  For  comparison,  the  structures  and  relative 
energies  of  the  stationary  points  for  the  5)^2  reaction  pathway, 
reported  previously  at  the  B97  1/ECP/d  level  of  theory,  are 
also  given  in  F^ure  2.  For  the  potential  energy  profiles  for  both 
the  proton  transfer  and  5^2  pathways,  the  first  stationary  point 
in  proceeding  fiom  the  reactants  to  products  is  the  HO”— 
HCH2I  hydrogen  bonded  prereaction  complex.  Also  shown  in 
Figure  2  is  the  energy  for  the  [CH3  I  OH]”  intermediate. 
With  zero  point  energies  included,  the  reaction  energies  for 
both  the  proton  transfer  and  pathways  are  in  quite  good 
agreement  with  eiqieriment.*®  The  B97  1/ECP/d  values  for 
these  proton  transfer  and  Sn2  reaction  energies  are  —4.2  and 
—64.9  kcal  mol”‘(with  zero  point  correction  energy  included), 
while  the  0  K  e:q)erimental  energies  of  reaction  are  —3.1  ±  1.8 
and  —66.4  kcal  mol”*.** 

As  shown  in  Rgure  2,  these  DFT  calculations  for  the  proton 
transfer  pathway  predict  the  existence  of  a  hydrogen  bonded 
prereaction  complex  HO”  HCH2I  (A)  and  a  transition  state 
[HO  H  CH2I]”  (TSl)  connecting  this  complex  with  the 
postreaction  complex  HOH  CH2I”  (C).  In  addition,  there  is  a 
second  transition  state  [ICH2  HOH]”  (TS2)  and  a  second 
postreaction  complex  CH2I”  H2O  (D).  For  the  prereaction 
complex  (a),  transition  state  (TSl),  and  postreaction  complex 
(C),  two  configurations  were  foimd  through  varying  the 
dihedral  angle  of  H— O— C— I.  If  we  use  cis  and  trans  to 
describe  the  structure  of  each  pair,  the  upper  set  of  energies  for 


A,  TSl,  and  C  has  H  and  I  in  the  trans  position,  denoted  as  A', 
TS',  and  C',  and  the  lower  set  has  H  and  I  in  a  ds  position, 
denoted  as  A",  TS",  and  C".  Although  the  stmctures  are 
(Merent  for  each  pair,  there  is  only  a  small  difiFerence  in  their 
energies.  The  DFT  energy  (without  ZPE  correction)  difier 
ences  are  0.02,  0.14,  and  1.22  kcal  mol”*  for  A,  TSl,  and  C, 
respectively. 

The  structures  for  the  A,  TSl,  and  C  stationary  points 
primarily  depend  on  the  distances  between  the  O  and  C  atoms 
and  a  H  atom  as  shown  in  Figure  3.  Table  5  gives  the  rQ^,  r^^ 
and  r^o  distances  for  the  stationary  points.  As  compared  to  the 
prereaction  complex  A,  the  corresponding  postreaction 
complex  C  has  a  shorter  OH  distance  and  longer  CH  and 
CO  distances.  The  second  transition  state  TS2  and  second 
postreaction  complex  are  found  by  rotating  the  H2O  molecule 
to  the  hydrogen  atom  side  of  the  CH2I  group,  as  compared  to 
the  structures  for  TSl  and  C.  Both  the  TS2  and  second 
postreaction  complex  structures  were  observed  in  the  trajectory 
simulations. 

C.  Probabilities  of  Atomistic  Mechanisms.  Reactions  1 
and  2  are  the  most  important  pathways,  and  in  previous  work** 
animations  were  used  to  identify  and  provide  detailed 
desorptions  of  their  atomistic  mechanisms  at  E„i  of  0.05,  0.5, 
1.0,  and  2.0  eV.  In  the  work  presented  here,  the  relative 
importance  of  the  different  mechanisms  versus  temperature  is 
considered. 

7.  Sfj2  Pathway.  As  observed  previously,**  two  direct 
mechanisms,  Le.,  rebound  and  strpping  and  several  indirect 
mechanisms  were  observed  for  the  5^2  pathway.  A  direct  fi'ont 
side  attack  mechanism  was  observed  in  the  previous  study  at  a 
collision  energy  of  2.0  eV,  but  this  mechanism  was  not 
observed  in  the  current  simulations  at  much  lower  collision 
energies.  The  indirect  mechanisms  observed  in  the  current 
study  include  the  roundabout,  formation  of  the  HO”  HCH2I 
hydrogen  bond  complex,  formation  of  the  CH3OH  I”  post 
reaction  compler^  central  barrier  recrossing  proton  exchange, 
and  couplings  of  these  events. 

Overall,  the  indirect  mechanisms  are  less  important  than  foe 
direct  mechanisms.  At  300,  400,  and  500  K,  the  firaction  of  foe 
indirect  mechanisms  is  0.47,  0.44,  and  0.41,  respectively.  The 
corresponding  re^ective  fractions  for  direct  rebound  and 
stripping  mechanisms  are  0.37,  0.34,  0.36,  and  0.16,  0.21,  0.23, 
dearly  showing  direct  rebound  to  be  more  important  than 
direct  stripping.  These  fi'actions  are  sununarized  in  Table  6. 

Figure  4  shows  foe  distribution  of  the  indirect  mechanisms. 
The  dominant  indirect  mechanism  involves  formation  of  foe 
HO”  HCH2I  complex.  Proton  exchange  is  the  second  most 
important  indirect  mechanism,  and  this  mechanism  is  always 
coupled  with  formation  of  foe  HO”  HCH2I  complex. 

A  quite  small  component  of  the  indirect  reaction  involves 
transient  formation  of  the  CH3OH  I”  postreaction  complex, 
as  seen  in  F^re  5.  The  number  of  such  tr^ectories  is  1  out  of 
311  at  300  K,  3  out  of  385  at  400  K,  and  2  out  of  301  at  500  K. 
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Proton  Transfer  Pathway 
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Figure  2.  Schematic  potential  profiles  for  the  OH”  +  CH3I  — ►  CHjF  +  H2O  protoi  transfer  and  OH”  +  CH3I  -»  CH3OH  + 1”  Sn2  reactions  at  the 
DFT/B97  1/ECP/d  level  of  theory.  The  energies  shown  ate  in  kcal  mol”*  and  are  relative  to  the  OH”  +  CH3I  reactants;  energies  listed  in  red 
specify  transition  state  structures.  Zero  pdnt  energies  are  not  included.  Two  configurations  were  found  for  the  first  transition  state  and  first  two 
intermediates  on  the  protrm  transfer  pathway  through  varying  the  dihedral  angle  of  H— O— C— L  If  we  use  cis  and  trans  to  describe  the  structure  of 
each  pair,  the  upper  set  of  energies  for  A',  TSl and  C'  has  H  and  I  in  the  trans  position,  and  the  lower  set,  denoted  as  A",  TS",  and  C",  has  H  and  I 
in  a  cts  positicxi. 


The  CH3OH  V  complex  is  a  stationary  point  on  die  Sn2 
pathway.  It  was  not  formed  in  the  previous  study*®  at  much 
higher  collision  energies  of  0.5,  1.0,  and  2.0  e'V.  Thougji  only  a 
very  small  percentage  of  the  current  trajectories  shows 
formation  of  the  transient  CH3OH  I”  complex,  this  study 
shows  that  at  lower  temperatures,  Le.,  collision  energy,  this 


complex  may  begin  to  participate  in  the  reaction  dynamics. 
This  complex  also  has  a  small  participation  in  the  OH”  +  CHjF 
Sn2  reaction.*® 

2.  CH2r  pathway.  The  proton  transfer  pathway,  i.e., 
formation  of  H2O  and  CH2l~  is  as  important  as  the  5^2 
pathway  at  300,  400,  and  500  K.  As  for  the  Sn2  pathway,  tiiis 
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Rgure  3.  Geometry  of  the  transition  state  [HO  H  CHjI]”  TSl  in 
Figure  2,  which  illustrates  the  distances  roH>  Tch,  and  r^o  h>r  Table  S 


Table  S.  Distances  Between  O,  H,  and  C  Atoms  for  the  Two 
Sets  of  Pre  reaction  Complexes,  Transition  States,  and  Post 
reaction  Complexes  for  the  OH“  +  CHjI  CHjI”  +  HjO 
Proton  Transfer  Reaction 


distance  (A)" 


v^per  set 

prereaction  complex 

A' 

'oh 

1.573 

'CH 

1.178 

'go 

2.728 

transition  state 

TSl' 

1.588 

1.174 

2.749 

postreaction  complex 

C' 

1.037 

1.771 

2.804 

lower  set 

prereaction  complex 

A' 

1AS3 

1.188 

2.728 

transition  state 

TSl' 

IA35 

1J03 

22532 

postreaction  complex 

C' 

1.033 

1.778 

2.811 

“The  r’s  are  defined  in  Figure  3. 


Table  6.  Simulation  Temperature  Dependence  of  the 
Fractions  of  the  Different  Atomistic  Reaction  Mechanisms 
for  the  OH-  +  CHjI  -r  CH3OH  +  I”  Sn2  and  OH’  +  CHjI 
— »  CHjI”  +  HjO  Proton  Transfer  Reactions 


temperature  (K) 

DR" 

DS" 

S,j2  reaction 

IND" 

300 

037  ±  0.04 

0.16  ±  0.05 

0.47  ±  OfM 

400 

035  ±  004 

031  ±  0.04 

0.44  ±  OfO 

500 

036  ±  0.05 

033  ±  0.05 
protnn-lransfer  reaction 

0.41  ±  0.04 

300 

0.06  ±  0.05 

033  ±  0.05 

0.71  ±  0D3 

400 

007  ±  004 

039  ±  0.05 

a64  ±  0.06 

500 

0.03  ±  0.03 

0.32  ±  0.05 

0.65  ±  0.03 

“DR,  direct  rebound;  DS,  direct  stripping;  IND,  indirect  mechanisms. 


proton  transfer  pathway  occurs  via  the  direct  rebound  and 
striping  mechanisms,  as  well  as  indirect  mechanisms.  Table  6 
displays  the  distribution  of  these  mechanisms.  The  indirect 
mechanisms  dominate,  and  the  fraction  ranges  from  0.64  to 
0.71.  Within  the  two  direct  mechanisms,  direct  striping  is 
more  important  than  direct  rebound. 

Several  types  of  indirect  mechanisms  were  assigned  for  die 
proton  transfer  pathway.  The  three  stationary  point  stmctures 
on  this  padiway,  observed  in  the  vast  m^ority  of  the  indirect 
trajectories  are  the  HO”  HCH2I-,  HOH  CH2I-,  and 
CH2r  H2O  complexes.  Identifying  HO”  HCH2I-  as  A, 
HOH  CH2I- as  C,  CH2I-  H2O  as  D,  roundabout  as  Ra,  the 
indirect  mechanisms  are  classified  as  Ra,  A,  C,  D,  and  the 
composites  of  the  complexes  A+C,  C+D,  A+D,  and  A+C+D. 
The  probabilities  of  these  eight  indirect  mechanisms  are  shown 
in  Figure  6.  In  at  least  90%  of  the  indirect  trajectories, 
structures  were  formed  for  HO”  HCH2I-  (A)  and/or 
042!“  H2O  (D).  It  is  noteworthy  that,  thou^  the  energy 
of  HOH  CH2r  (C)  is  lower  than  that  of  CH2r  H2O  (D), 


Mechanism  Type 

Figure  4.  Fractions  of  the  individual  atomistic  reaction  mechanisms  for 
the  OH”  +  CHjl  -♦  CHjOH  +  F  Sn2  reaction  at  diSerent 
temperatures:  DR,  direct  rebound;  DS,  direct  stripping;  A,  HO” 
HCH2I  prereaction  complex;  hr,  barrier  recrossing;  PE,  proton 
exchange;  and  B,  CH3OH  F  postreaction  complex. 
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Figure  S.  Frames  of  a  sample  CH3OH  F  postreaction  complex 
fiirmation  tiajectcay  leading  to  the  CH3OH  +  F  S|.j2  products  via  an 
indirect  mechanism. 


Mechanism  Type 

Figure  6.  Fractions  of  the  individual  atomistic  reaction  mechanisms  for 
the  OH"  +  CH3I  -*  CH2F  +  H2O  proton  transfer  reaction,  at 
different  temperatures:  DR,  direct  rebound;  DS,  direct  stripping;  Ra, 
roundabout;  A,  HO”  HCHjI  prereaction  complex;  C,  HOH  CHjF 
postreaction  complex;  and  D,  GHjI*  H2O  postreaction  complex. 


D  is  favored  before  H2O  and  CH2F  separate.  This  is  because 
after  their  formation  both  CH2F  and  H2O  have  adequate 
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energy  to  rotate  around  each  other,  giving  rise  to  a  long  lifetime 
for  the  CH2I”  H2O  complex. 

Several  other  relatively  unimportant  mechanisms  occurred 
for  a  very  few  of  the  trajectories,  which  are  coupled  with  the 
above  mechanisms.  For  one  of  these  mechanisms,  CH3OH  and 
1“  were  first  formed,  then  associated  to  form  CH3OH  T,  and 
through  a  rearrangement  and  proton  transfer  formed  CH2l~  + 
H2O.  For  another,  proton  transfer  occurred  twice,  and  die 
successful  transfer  forming  products  was  completed  by  an  H 
atom  fi'om  CH3I  different  than  the  one  for  the  first  proton 
transfer.  The  remaining  type  is  OH“  and  CH3I  collision  to  form 
[CH3  I  OH]~,  which  existed  for  a  while  before  dissociating  to 
CH2r  and  H2O. 

D.  Impact  of  OH~  Rotational  Excitation.  Simulations 
were  also  performed  to  study  die  role  of  OH“  rotational 
excitation  on  the  reaction  dynamics.  The  previous  simulations 
at  a  collision  energy  of  O.OS  eV,  with  fixed  OH”  rotational  and 
vibrational  quantum  numbers  /  =  3  and  n  =  0,  and  CH3I 
rotational  and  vibrational  temperatures  T,  =  330  K  and  T,  = 
130  K  were  repeated,  but  with  /  equal  to  0  and  6.  A  total  of 324, 
75S,  and  325  trajectories  were  calculated  for  /  of  0,  3,  and  6, 
respectively.  The  reaction  probability,  anion  product  ratios,  and 
mechanisms  for  the  S142  pathway  were  analyzed  as  a  function  of 
/• 

As  shown  in  Table  7,  the  overall  reaction  probability  and  the 
probabilities  for  the  different  pathways  have  at  most  only  a 


Table  7.  Simulation  Probabilities  for  the  OH”  +  CH3I 
Reaction  Pathways  versus  OH”  Rotational  Quantum 
Number 


7(OH-) 

0 

3 

6 


probability 


total  reaction 
0.58  ±  0.03 
0.49  ±  0.02 
0.56  ±  0.03 


CH3OH  +  r 
0.33  ±  ao3 
0.28  ±  0.02 
0J3  ±  a03 


CHjr  +  HjO 
0.22  ±  0.02 
0J8  ±  0.01 
020  ±  0.02 


[CHj-I-OH]- 
0.03  ±  aoi 
0.03  ±  0.01 
0.03  ±  0.01 


small  dependence  on  /.  There  is  an  indication  that  the  overall 
reaction  probability  is  sl^hdy  smaller  at  /  =  3,  but  this  is  a 
minor  effect.  Consequendy,  the  anion  product  ratio  is  similar 
for  the  three  /  values.  As  shown  in  Table  8,  the  ratio  for  I”: 


Table  8.  Simulation  Anion  Product  Ratios  versus  the  OH' 
Rotational  Quantum  Number 


J  (OH-)  CHiJ- 

0  028  ±  0.04 

3  027  ±  0.06 

6  a35  ±  0.04 


[CHj-l--OH]-  r 

0.06  ±  0.01  026  ±  0.05 

0.06  ±  0.01  027  ±  0.06 

0.06  ±  0.01  029  ±  0i)5 


striping  mechanisms  are  dominant.  However,  as  /  is  increased, 
the  fiaction  of  the  direct  mechanisms  increases  and  the 
corresponding  fiaction  of  the  indirect  mechanisms  decreases  as 
0.43  ±  0.06,  0.39  ±  O.OS,  and  0.36  ±  0.06  for  /  of  0,  3,  and  6, 
respectively.  The  fi-actions  for  the  indirect  and  direct 
mechanisms  are  summarized  in  Table  9. 


Table  9.  Simulation  Fractions  of  the  Different  Atomistic 
Reaction  Mechanisms  for  the  OH”  +  CH3I  -+  CH3OH  +  I” 
S„j2  Reactions  versus  OH”  Rotational  Quantum  Number 


7  (OH”)  DR  DS 

0  029  ±  004  027  ±  004 

3  023  ±  0.06  028  ±  0.05 

6  029  ±  006  025  ±  004 


IND 

0.43  ±  0.06 
029  ±  0.05 
026  ±  0.06 


“DR,  direct  rebound;  DS,  direct  stripping  IND,  indirect  mechanisms. 


Although  the  overall  trend  is  for  the  direct  Sm2  mechanisms 
to  increase  with  increase  in  /,  this  is  not  the  case  for  the 
individual  direct  mechanisms.  The  probability  of  the  direct 
rebound  mechanism  increases  with  increasing  /,  while  die 
probability  of  the  direct  striping  mechanism  remains  the  same 
with  change  in  /.  The  increased  fraction  of  the  direct 
mechanisms  with  increase  in  /  is  contributed  to  the  direct 
rebound  mechanism. 

The  probabilities  of  the  different  indirect  S2j2  mechanisms 
are  not  affected  by  /.  As  shown  in  F^re  7,  formation  of  the 


Mechanism  Type 

Figure  7.  Fractions  of  the  individual  atomistic  reaction  mechanisms  for 
the  OH”  +  CH3I  -+  CH3OH  +  P  5^2  reaction  versus  the  OH” 
rotational  quantum  number:  DR,  direct  rebound;  DS,  direct  stripping; 
Ra,  roundabout;  A  HO”  HCH2I  prereaction  complex;  br,  barrier 
recrossing;  PE,  proton  exchange;  and  B,  CHjOH  I”  postreaction 
complex. 


CH2r:  [CH3  I  OH]”  is  0.56:0.38:0.06,  0.57:0.37:0.06,  and 
0.59:0.35:0.06  for  /  of  0,  3,  and  6,  re^ectively.  With  [CH3  I 
OH]”  not  included  as  a  product  (see  IVA.1  Simulation 
Results),  the  respective  rate  constants  are  2.2  ±  0.2  X  10”’,  1.9 
±  0.1  X  10”’,  and  2.1  ±  0.1  X  10”’  cm^s”'  for  J  ofO,  3,  and  6, 
respectively.  In  summary,  these  simulations  indicate  that  change 
of  the  OH”  rotational  quantum  number  does  not  have  a 
significant  effect  on  the  probability  of  the  overall  reaction,  nor 
for  the  different  pathways. 

A  detailed  analysis  of  the  dynamics  for  the  S;,j2  pathway 
shows  that  the  value  of  the  OH”  rotational  quantum  number 
does  influence  the  atomic  level  mechanisms.  As  for  the  results 
presented  above  versus  temperature,  the  direct  rebound  and 


HO”  HCH2I  complex  is  the  dominant  mechanism.  The 
proton  transfer  mechanism  is  second,  followed  by  barrier 
recrossing,  postreaction  complex  formation,  and  roundabout  as 
lesser  important  mechanisms.  In  sutmnary,  the  value  of  /  affects 
the  overall  distribution  of  indirect  and  direct  mechanisms  but 
not  the  individual  indirect  mechanisms. 

V.  CONCLUSION 

Both  experiments  and  simulations  find  that  the  OH”  t  CH3I 
reaction  for  temperatures  between  300  and  500  K  occur  via  a 
Sn2  pathway  forming  CH3OH  +  I”  and  a  proton  transfer 
pathway  forming  CH2I”  +  H2O.  The  relative  importance  of 
these  two  pathways  are  similar;  the  experimental  F  fiaction  of 
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the  anion  yield  decreases  from  0.64  to  0.46  for  a  temperature 
change  from  210  to  500  K,  vdiik  the  simulation  fractional  I~ 
yield  decreases  as  0.56  ±  0.06,  0.55  ±  0.04,  and  0.50  ±  0.05  for 
the  respective  temperatures  of  300,  400,  and  500  K.  The 
reaction  intermediate  [CHj  I  OH]“  is  present  in  the 
simulations  when  the  trajectories  are  ended  at  2.5  ps,  but  this 
^edes  dissociates  primarily  to  the  reactants  CH3I  +  OH“  at 
longer  times  and  is  not  observed  in  the  experiments. 

The  atomistic  mechanisms  for  the  OH“  +  CH3I  reaction  are 
at  most  only  weakly  affected  by  the  reactant  temperature.  For 
the  Sn2  pathway,  the  direct  mechanisms  are  sli^tly  more 
important  than  the  indirect  mechanisms,  and  reboimd  is  the 
dominant  direct  mechanism.  The  indirect  mechanisms  are 
more  important  than  direct  mechanisms  for  the  proton  transfer 
pathway,  and  stripping  is  the  dominant  direct  mechanism. 

Rotational  exdtadon  of  OH“  does  not  affect  the  overall  OH“ 
+  CH3I  reaction  probability  nor  tfie  branching  between  the  Sn2 
and  proton  transfer  pathways.  However,  it  may  affect  the 
probabilities  and  dynamics  of  the  atomistic  mechanisms  for 
these  pathways. 
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